Spin torque and waviness in magnetic multilayers: a bridge between Valet-Fert theory 

and quantum approaches 
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We develop a simple theoretical framework for transport in magnetic multilayers, based on 
Landauer-Buttiker scattering formalism and Random Matrix Theory. A simple transformation 
allows one to go from the scattering point of view to theories expressed in terms of local currents 
and electrochemical potential. In particular, our theory can be mapped onto the well established 
classical Valet Fert theory for collinear systems. For non collinear systems, in the absence of spin-flip 
scattering, our theory can be mapped onto the generalized circuit theory. We apply our theory to 
the angular dependance of spin accumulation and spin torque in non-collinear spin valves. 

PACS numbers: 72.25.-b; 72.25.Ba; 75.47.-m; 85.75.-d 
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The discovery of Giant Magnetoresistance (GMR) [l[ 
and its subsequent counterpart, the current induced spin 
torque 0] , was at the origin of a new field which aims at 
controlling the magnetization dynamics of small metal- 
lic devices through standard electronics. The theory of 
transport in those systems is now well developed and in- 
cludes a number of approaches that range from classical 
Valet Fert (VF) theory in the diffusive regime [3J , Boltz- 
mann equation Q| to quantum approaches (original [f| 
and generalized |6j) circuit theory, random matrix theory 
(RMT) for the scattering matrix [7| , and ab-initio based 
models [3, 0] . Many connections exist between these dif- 
ferent approaches. One popular route flo| starts from 
the Keldysh Green's function formalism. In the quasi- 
classical approximation it yields the Boltzmann equa- 
tion and the (VF) diffusive equation. A different strat- 
egy where local equilibrium is assumed only in certain 
points leads to circuit theory. An alternative route is the 
Landauer-Buttiker formalism which expresses the prob- 
lem of transport inside a quantum conductor as a scat- 
tering problem where electrons undergo various trans- 
mission/reflection events as they go through the system. 
This approach is well suited for coherent system and is 
equivalent to the Keldysh approach. However, the clas- 
sical concepts of chemical potential or local equilibrium 
do not arise naturally in the scattering approach so that 
classical intuitions do not easily transfer into its language. 

In this letter we take the scattering formalism as our 
starting point and develop a theory which fully captures 
VF and (generalized) circuit theory. Our theory (here 
after referred as C-RMT for Continuous Random Matrix 
Theory) can be tabulated by the same set of (experimen- 
tally accessible) parameters as VF [ll|. On the other 
hand, it properly includes Sharvin resistance and allows 
for non collinear and even a one-dimensional texture of 
magnetization (i.e. domain walls). We apply C-RMT to 
the discussion of the angular dependance of spin torque 
("waviness" [HI, El) in asymmetric spin valves. 
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FIG. 1: (a) Schematic of a spin valve with two ferromagnetic 
layer Fa and Fb whose magnetization makes an angle 9. (b) 
and (c): geometric construction of the spin torque, see text. 
J a, Jn and Jb are the spin currents along the valve while ta 
and tb are the torque on the two magnetic layers. 



Two quantum conductors in series. Before introducing 
C-RMT, we start with an instructive example of trans- 
port in a non-magnetic metal. C-RMT is a mere exten- 
sion of the concepts developed below to properly take into 
account the electronic spins. Classically, the electronic 
transport in a system A is described by its resistance TZa- 
When two systems A and B are put in series, the total 
resistance is described by Ohm law TZab = T^-A + TIb- 
In a coherent quantum system, the transport properties 
can be described by the amplitudes for reflection and 
transmission, and the conductor is characterized by its 
scattering matrix Sa, 
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where ta and tA it' a an d r' A ) describe the reflec- 
tion/transmission amplitudes to the right (left) of the 
sample. When putting two such conductors in series, one 
has to add the amplitudes for all possible processes. For 
instance t A B = tBtA+tB(rAr' B )tA+tB(r A r' B ){rAr' B )tA + 
. . . includes the direct process where the electron is trans- 
mitted by A and then by B {t B tA), the process where the 
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electron is transmitted by A, then reflected by B and A, 
and then transmitted by B (tBTAr' B tji) and so on. The 
previous geometrical series can be resumed into, 
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The conductance of such a system is given by the Lan- 
dauer formula qab = (e 2 /h)TAB where Tab = \tAB\ is 
the probability for an electron to be transmitted. Let 
us now suppose that many channels N c h of propaga- 
tion are present in the system, so that the transmission 
probability has to be averaged over those. The product 
r' B rA = \r' B rA\z contains a phase z = e 1 ^ which is al- 
most random from channel to channel, so that with very 
good approximation the total transmission is given by its 
average over z. Straightforward integration gives, 



Tab =T b (1-R a R' b )- 1 Ta 



(4) 



where Ra = \ta\ 2 — 1 ~ Ta is the reflection probability. 
Eq. flU) can be recast into 1 /T AB = 1 /T A + 1 /T B - 1 which 
leads to 1Z ab — 71 a + 71b — 7t s h where the Sharvin resis- 
tance is given by 7Z s h = h/ (e 2 7V c h). The above equations 
call for a number of comments, (i) Addition law for resis- 
tances is very close to Ohm law except for the presence of 
the Sharvin (or contact) resistance, but for classical con- 
ductors, TZa(b) ^ 7Z s h one recovers Ohm law. (ii) Eq.Q 
has a structure completely similar to Eq.Q. In fact, it 
can be obtained directly in a way similar to Eq.([2|) by 
adding not the amplitudes but the probabilities for the 
various processes, (iv) We note that Ohm law is derived 
here using purely quantum mechanical concepts for fully 
coherent conductors. In particular there is no need for 
a well defined chemical potential (local equilibrium) in 
between the two conductors A and B even though every- 
thing happens as if such a chemical potential was defined. 

RMT. We now proceed with the extension of the above 
ideas to magnetic systems where the S matrix now in- 
cludes a spin grading. Our starting point is an extension 
of RMT that was introduced in Ref 0| ■ The theory has 
a structure almost identical to the scattering approach 
except that a conductor (or a part of it) is now described 
by 4 x 4 "hat" matrices f,f',t and t' defined in term of 
the (4A ch x 4A ch ) S matrix, 
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with similar definitions for f, i and V . The main result 
of Ref. p]] is that the addition law for the "hat" matrices 
is exactly given by Eq.((2]) and (J3j) except that now the 
transmission and reflection matrices are to be replaced 
by their "hat" counterparts. The conductance is given 



by .9 = (hi + hi + hi + hi) /Tlsh = (1/^sh) J2aa> T <?<?'i 

where T aa > = (1 / 'N<&) r Tit rTrT iV cr(J , is the probability for an 
electron with spin tr' to be transmitted with a spin a. 
The "hat" matrices for the interfaces between two metals 
(say copper and cobalt) can be obtained from ab-initio 
calculation or tabulated from the experiments 

(see below). Once the "hat" matrices are known for one 
direction of the magnetization, they can be rotated in 
arbitrary direction using rotation matrices [7| . 

C-RMT. To obtain the "hat" matrices for the bulk 
parts of the metallic layers, we need to proceed further. 
Assuming we know the "hat" matrix S(L) for a system 
of size L, we can add an infinitely small layer of size 8L 
and compute S(L + 5L) from the knowledge of S(5L) and 
the addition law for "hat" matrices. The generic form for 
S(SL) is given by, 



t = 1 - A*5L 



A r 5L 



(0) 



where the two matrices A* and A r entirely characterize 
the bulk properties of the material. Expanding the addi- 
tion law Eq.([2|) in 5L provides, 

dr/dL = A r - A*f - ?A* + fA r f (7) 
di/dL = -KH+rPSt (8) 

The various elements entering in the " hat" matrices (say 
t) are vastly inequivalcnt. The main elements are the 
probabilities T aa i . Secondly comes the element on the 
diagonal, the so called (complex) mixed transmission 
T mx = (l/JVoi,)Trt-["[-i| i which measures how a spin trans- 
verse to the magnetic layer can be transmitted through 
the system. These elements are usually believed to be 
rather small TH, 3] in magnetic systems, but can play 
a role in non-collinear configurations nevertheless. Last, 
the other elements involve some coherence between spin- 
flip and non spin-flip processes and are likely to be even 
smaller. In the basis parallel to a layer magnetization, 
they can be disregarded. We parametrize a bulk layer by 
four parameters T^, T| T s f 

( r T + r sf 
o I 
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(9) 



l a = 1/Y a is the mean free path for spin a. In a ferro- 
magnet, L mx = l/l± + I/Il where l± is the penetration 
length of transverse spin current inside the magnet while 
II is the Larmor precession length. Those lengths, which 
are roughly equal, are the smallest characteristic lengths 
with typical values smaller than lnm. In a normal metal, 
L T = r x = r and T mx = T + 2r sf so that the "hat" matri- 
ces remain invariant upon arbitrary rotation of the spin 
quantization axis. A r is given by the same parametriza- 
tion as Eq.Q with T s f being replaced by — T s { (in order 
to fulfill current conservation) and neglecting r mx for fer- 
romagnets (as the mixing conductance is essentially of 
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ballistic origin). This completes the formulation of the 
theory. Eq.© and Eq.® can be integrated to obtain 
the "hat" matrices of the bulk parts. A given multilayer 
is then constructed by using the addition law Eq. © and 
([3]) for the various bulk layers and the corresponding in- 
terfaces. 

Link with Valet Fert theory. Let us introduce the 4- 
vector P±0) = {P±,h p ± ,mx) -P±,mx> P±,l where P + ^{x) 
(P-,l) is the probability to find a left (right) moving elec- 
tron with spin up (down) at point x. The addition law 
Eq.© and (O (for "hat" matrices) is equivalent to state 
that P±(xi) and P±(iC2) on two sides of a conductor A 
are related through its "hat" scattering matrix Sa as, 



P-(*i) 

P+(*2) 



= S A 



P+(*i) 

P-(»2) 



(10) 



In the physical picture where one represents the scatter- 
ing processes as random events with certain transmission 
or reflection probabilities, the above equation has a sim- 
ple interpretation when one focuses on its first and fourth 
raw -P±,t and P±,i- it accounts for the conservation of 
probability in the scattering events, i.e. it is the Master 
equation of the underlying Brownian motion undertaken 
by the incident electrons. Let us now introduce two new 
4- vectors, and fi(x) defined as, 

j(aO = [P + (z)-P_(z)]/(e7e 8h ) (11) 
l*(x) = [P+(x) + P-(x)}/2 (12) 

(with e < 0). Now, using the parametrization Eq.® 
and Eq.© for infinitely thin layer, and writing Eq. (|10p 
in term of j(x) and n(x), we arrive at, 

ja = -l/(eT a TZ sh ) d x pb a (13) 
d x j a = 4L sf /(e^ sh ) - fi a ] (14) 

which are precisely the VF equations Hence, for a 
collinear system, C-RMT simply reduces to VF theory. 
In its original form however, VF theory does not account 
for the presence of Sharvin resistance. Here the bound- 
ary conditions are given by the Landauer formula: the 
presence of a potential drop eV between the reservoirs 
located at, say, x = and x = L imposes P +a (0) = eV 
and P^ a (L) = which translates into, 

/v(0) + (e7W2) j a (0) = eV (15) 
(i a (L) - (eU sh /2) j a {L) = (16) 

These mixed boundary conditions allows VF theory to 
properly include the Sharvin resistance of the system and 
correspond to adding half of the Sharvin resistance on 
both extremities of the sample. 

Link with generalized circuit theory. Let us now con- 
sider Eq. (fT0|) for a conductor whose transmission and re- 
flection matrices are purely diagonal, i.e. without any 
spin-flip scattering. Let us further suppose that i? mx 



^0 
-0.1 



_«/ ^"1000 


COo 






Dooooopooood 


i 



Cu 



10 



1000 x(nm ) 1020 



Cu 5 Au m 




02 0.4 0.6 0.8 1 

9/k 



15 



FIG. 2: a) Spin accumulation in the middle of a spin valve 
A — CuioooCosCuioPysCu^Au^oo (thickness in nm) for dif- 
ferent angle 9 = (circles), 7r/4, tt/2, 3n/4 and tt (squares). 
Symbols stand for VF calculations while lines correspond to 
C-RMT. b) torque tb(0) (per total current j) on the Py-layer 
of A for various Py-thickness Lp y from 0.5nm (thick) to 15nm 
(dashed) . c) Stability angle 8* as a function of Lp y for A (cir- 
cles), B = CuioooCosCuioPyLp y Cuiooo with 8 = (squares) 
and B with ScoCu = &CoPy = 0.25 (diamonds). 



might be non zero but T„ 
form 

3<i{xi) = j<r(x 2 ) 



0. Then, Ea.([TD]) takes the 
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T, 



e1Z s h 1 - T, 



?—[(i a (xi) - fj, cr (x 2 )^17) 



9 1 — Z? 1 * 2 

Jmx(Sl, 2 ) = ± 



where T a = T\ a + T^ a is the total transmission of an 
incident electron with spin a and R]£ are the mixing 
reflections from left to left (i? 1 ^) and right to right 
(i?mx)- Eq. flTTf and Eq. (fl~8|) define the generalized cir- 
cuit theory [6j so that in the absence of spin-flip scatter- 
ing C-RMT and generalized circuit theory are completely 
equivalent. In fact, the renormalization coefficients of 
generalized circuit theory Q were chosen such that the 
calculation of the conductance with RMT and general- 
ized circuit theory fully agree with each other. We find 
that the point of view of scattering taken in this letter 
is fully equivalent to the alternative view in term of lo- 
cal current and chemical potential (VF, circuit theory), 
and one can change from one to the other simply using 
Eq. (fTTj) and Eq. (fT2"|) . We note that this mapping is quite 
general and can be extended to include for instance in- 
terface spin-flip scattering or superconductivity [13] ■ 

Tabulation of CRMT with Valet Fert set of parame- 
ters. As C-RMT and VF are equivalent, we can use the 
huge corpus of experimental data that has been inter- 
preted within VF to parametrize C-RMT. The Valet- 
Fert resistivities p-m) = 2p*(l =F (3) and spin-diffusion 
length Z s f are in a one-to-one correspondence with C- 

rmt: i/i sf = 2VWr T + r x , p = {Y i - r T )/(r T + ro 

and p*/TZ sh = (r T + r x )/4. Similarly, within VF, inter- 
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faces are described by a set of three parameters, r^,^ — 
2r b *(l =F 7) and 5 which correspond to an effective layer 
of thickness d which is taken to be infinitely thin while 
keeping the parameters S = d/l s f and r a — p a d constant. 
Figf2](a) shows the spin accumulation profile of a Co/Py 
spin valve calculated with C-RMT. For the parallel and 
anti-parallel configuration we also show an (independent) 
Valet-Fert calculation using the CNRS-Thales software 
developed by one of us (HJ). Both match perfectly as 
expected. 

Application to spin torque: Wavy or not wavy? We 
now turn to an application of C-RMT to the calculation 
of spin torque in realistic spin valves, taking int account 
spin-flip scattering both in the bulk and at the inter- 
faces. We consider the valve depicted in FigQ] (a) and 
note J4, J/v and Jb the spin currents just before, in be- 
tween and after the two magnetic layers Fa and Fb ■ For 
non collinear magnetization, spin current is not conserved 
and the spin torque on Fa and Fg is defined as, 



ta = J A ~ Jn, tb — Jn — Jb 



(19) 



Let us start with a simple geometric construction that 
allows to get a physical picture for the torque in a rather 
general way. To do so, we need two hypothesis: (i) the 
mixing transmission are small (it is the case for metallic 
magnetic layers) so that J a and Jb are parallel to the 
magnetization of Fa and Fb respectively, (ii) the sys- 
tem is thin enough for spin-flip scattering to be ignored 
in the active region so that ta (tb) is perpendicular to 
J A {Jb)- The construction goes as follows, see FigfT] (b) 
and (c): first we plot J a and Jb which make the same 
angle 9 as the magnetization of the respective magnetic 
layers. Then, we note that r tot = ta + tb = J a — Jb 
does not depend on the unknown Jn and points from 
the tip of Jb to the tip of J a- The construction of 
the torque is then straightforward: the two vectors ta 
and tb are chosen such that they are perpendicular to 
their respective layer and their sum goes from the tip 
of Jb to the tip of J a- This simple construction gives, 
in particular, the sign of the torque as a function of the 
angle 9. We find that when, say, Jb > J a the torque 
on the layer with the highest polarization (Fb) can be- 



come wavy [12|, |13|, i.e. instead of favoring the parallel 



or antiparallel configurations, the torque stabilizes (or 
destabilizes depending of the direction of the current) a 
configuration with a finite angle 9*. The critical angle 
9* where the torque vanishes verifies \Ja/Jb\ — cos 9*. 
On the other hand, at small angle, one has the follow- 
ing developpement \Ja/Jb\{0) = 1 — r/9 2 /2 + ... (Current 
conservation imposes J a = Jb at 9 = and the ratio is 
an even function of 9) so that "waviness" is found when 
77 > 1. As rj = for a symmetric structure, it means that 
a finite asymmetry is needed to enforce waviness. 

Without spin-flip scattering, we find rj — Yibt*b ~ 
1AT* A + (iB - 7A)? , A r s/'^sh]/(7sfB + lAr* A ) where the 



effective parameters ja and r* A include both the inter- 
face and bulk properties of layer A. More generally, the 
crossover between normal and wavy can be discussed by 
looking at the small angle expression derived by Fert et 
al (Eq.(5) in |l8[, we omit the ballistic corrections). It 
can be obtained by relating the spin accumulation in the 
spacer for 9 -C 1 to spin current and spin accumulation 
at 9 = 0, and then applying Eq. ^Tg]) . It reads, 



dT B 

~~dT 



0=0 



— - — \e=o + n ^ — |e=o 



2eK, 



sli 



(20) 



When, for instance, one crosses the assymetry border 
from J a > Jb > to Jb > J a > 0, the spin accumu- 
lation, proportional to the gradient of the spin current, 
changes from negative to positive and the second term 
in Eq. (|2"0|) becomes neagative and begins to compensate 
the first one. The crossover from normal to wavy occurs 
when, by a further increase of assymetry, the spin accu- 
mulation term wins and reverse the sign of dTs/d9\Q = Q 
(see Gmitra and Barnas for an extensive discussion 
of the normal to wavy crossover). 

Typical examples of our numerical results are pre- 
sented in Figf2]for Co/Cu/Py samples in which the asym- 
metry comes from the short l s { and large polarization and 
resistivity of Py. Starting from a small value of Lp y , an 
increase of the assymetry and finally a crossover to wavy 
(9* =^ 0) can be obtained by increasing Lp y as shown in 
Figf2] (b) and (c). By comparing the curves in samples 
with and without Au on the right of the valve, one sees 
that the short l s { of Au in a layer close to Py tends to 
increase the asymmetry and the waviness. On the other 
hand, interface spin-flip is found to favour a normal spin 
torque. In the experimental results of Boulle et al [13j 
a wavy behaviour was found for Cu/Co/Cu/Py/Cu/Au 
structures with equal thicknesses (8 nm) for Co and Py. 
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ful discussions, S. Petit- Watelot for his insights on the 
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